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COMPUTER SIMULATION 

FLIGHT ACCELERATION FACILITY ENTRY MONITORING SYSTEM STUDY 


INTRODUCTION 

The Engineering and Development Directorate was asked to support 
the Flight Crew Operations Directorate in a closed loop Apollo entry 
simulation using the Flight Acceleration Facility (centrifuge). The 
objective of the simulation was to provide the flight crew a real time 
simulation to use the entry monitor system for lunar return entry 
conditions while subjecting the crew to realistic accelerations. At 
the time of the request, the flight acceleration facility of the Crew 
Systems Division did not he.ve the computing capacity for a simulation of 
this size and the Guidance and Control Division was asked to aid in the 
work. The result was a Joint effort of the Crew Systems Division and the 
Guidance and Control Division with the Crew Systems Division having overall' 
responsibility for the job as well as supplying the centrifuge and the 
normal complement of supporting personnel and equipment. The Guidance 
and Control Division supplied the equations of motion, two analog computers, 
and programed and operated the analog computers. The size of the task 
was such that the third computer had to be borrowed from the Computation 
and Analysis Division. Within the Guidance and Control Division, the 
Control Requirements Branch developed the vehicle equations of motion and 
aided in the checkout of the simulation. The Hybrid Computation and 
Simulation Branch programed the equations, checked the simulation, and 
maintained the computers. 

The purpose of this paper is to describe the analog simulation of * 
the EMS equations of motion and how it interfaced with the overall system. 
The operation of the centrifuge and the programing of the centrifuge 
control computers will not be discussed in any detail. 




SYMBOLS 


Applied acceleration along 3^, Y^, and Z^, g’s 

t 2 

Applied acceleration along X^, Y^, and Z^, ft/ sec 

Aerodynamic axial and normal force coefficients and 
pitching moment coefficient, respectively 

Slopes of aerodynamic C A , (? N , and C m versus cosine r\ 
curve, respectively 

0 

Spacecraft reference length, ft 

2 

Acceleration due to gravity, ft/sec 

Spacecraft altitude above sea level and altitude 
rate, ft, and ft/sec., respectively 

Moment of inertia about 

2 

respectively, slug-ft 
Gravitational gain factor 

Ordinate value of C A , C. T , C versus, cosine n curve, 

A U m 

respectively. 

RCS jet rolling, pitching, and yawing moments, 
respectively, ft -lb. 

Spacecraft mass, slugs 

Roll, pitch, and yaw rates and accelerations about 
3^, Y^, and Z^, respectively, deg/sec and deg/sec^ 

g 

Dynamic pressure, lb/ft 

Reference radius of the earth, ft 

x . 




and 




axes 


3 


U t V, V 

• • • 
U, V, w 


V, 


T 


eg 


tano< 
tano< A 
Cosine tj 

e.f, 0 
4,*, 0 


e 

k 

K 


2 

Spacecraft reference area, ft 

Body components of inertial velocity, ft/sec 

Components of body acceleration, ft/sec^ 

Total velocity, ft/sec 

Distance from spacecraft center line to c.g. along 

Tangent of trim angle of attack (SCS) 

Tangent of trim angle of attack (aero) 

Cosine of total angle of attack 

Fuler angles and Euler angle rates, respectively, 
deg 8nd deg/sec 

Atmospheric density, slug/ft 3 

I 

Aerodynamic resolution angle, deg 

Spacecraft roll angle about the velocity vector, deg 
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DESCRIPTION OF SIMULATION 
General 

The computers used for this simulation were three AD80 analog 
computers and two DDP24 digital computers. The three analog computers 
were used to simulate in real time the vehicle dynamics during the 
entry phase of a lunar return trajectory. The analog outputs, equation 
solutions, were then fed via A-D converters to the DDP24 digital control 
computer and to the DDP24 digital aerodynamic computer. The aerodynamic 
computer took the analog information and calculated the required angular 
velocity for the centrifuge arm and gondola attitude position. This 
calculated data was then transmitted Via the DDP24 control computer to 
the centrifuge servo control electronics which drove the centrifuge aim 
motor end gondola gimbal to properly simulate the acceleration forces 
acting on the command module. Also, the vehicle attitude signals (from 
the analog) were fed to the cockpit instruments via the DDP24 control 
computer. To close the loop, the pilot, using the attitude information 
end EMS (entry monitor system) information in addition to other informa- 
tion, sent commanded roll signals via the RHC (rotational hand controller) 
to the analog computer. Figure 1 is a flow diagram of the total system. 

All computers, analog and digital, were in a first floor laboratory 
end the centrifuge and centrifuge electronics in a separate room in the 
Flight Acceleration Facility, Building 29. The two DDP24 digital computers 
were programed and operated by System Test Branch personnel, Crew Systems 
Division. Besides controlling the centrifuge and cockpit instrumentation, 
the DDP24 computer had overall control of the total simulation node status, 
that is, the analog modes (initial condition, hold, and operate' were 
controlled by signals from the DDP24 computer, This will be discussed 
in more detail in "special circuits." For more detail on the operation 
of the centrifuge and digital computers, see reference 1. 


Description of Vehicle 


Vehicle mass properties .- The command module mass, inertias, and 
center of gravity were all held constant. The vehicle mass characteristics 
were as follows: 


Mass 

I 


III 
xy' yz> zx 


= 402.6 slugs 
- 6250.1 slug -ft 2 

= 5851.2 slug-ft 2 

«= 5851.2 slug-ft 2 

= 0.0 
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The coordinates of the vehicle center of gravity expressed in terms 
of the standard Apollo coordinate system were X og - 104l.6 inches, 

Y * 0.0 inches, and Z *5-57 inches. This center of gravity gave 
eg eg 

the vehicle a llft-to-dreg ratio of .28 which vac nominal for the Apollo 
101 command module. 


Vehicle RCS thrust characteristics .- During entry from .05 g altitude, 
the only control thrust 3va liable is from the CM RCS jets. For this simu- 
lation the rotational moments resulting from CM RCS jets were as follows: 


* i j ■ «- im 

♦“j ■ 513 «-“■ 

+Nj ■ 478 ft -lbs 


-Lj - -452 ft-lbs 
-Mj --340 ft-lba 
-Nj ■= -480 It -lbs 


Rotational moments due to jet interaction were not considered. 

Vehicle control system .- The command module reaction control system 
used to command the rotational moments consisted only of rate -command 
control mode with a rate deadband of +2 degrees/sec and a hysteresis 
of .2 deg/sec. Inputs into the control system were the pilot's commanded 
rotational rates via the rotational hand controller. The commanded rates 
were proportional to RHC deflection with a maximum command of 20 deg/sec. 

A block diagram of the control system is shown in figure 2. 

Vehicle aerodynamic characteristics .- The command module which 
was simulated to have the aforementioned mass properties and an LOD 
(lift-to-drag ratio) of .280 was aerodynamics lly stable for earth entries 
at a trim angle of attack of 18.7 degrees. Angle of attack (tj) is defined 
as the acute angle between the total velocity vector and the («X) body 
axis. The equation defining angle of attack is 

T) * arc cosine (~) where: V^ = total velocity 

U = projection of V r ^ on X body axis 
Figure 3 is a diagram showing the angle of attack, tj. 

The axial and normal force coefficient slopes and pitching moment 
coefficient slope were considered linear and a function of cosine tj and 
not as a function of Mach number. The equations used to define the 
aerodynamic coefficients are the following: 



■ C A 

Cos 

+ K 

A 

n 


= C N 

Cos_ 

T 

+ K 

it 

o 

s 

1 

Cos^ 

T1 

+ K 


2 

3 

k 
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where: Ordinate value C A vs. Cos r\ curve, K 2 * 1.05 

Ordinate v«lue vs. Cos r\ curve, * ,112k 

Ordinate value C vs. Cos n curve, Kl * .3429 

ro h* 

Axial force coefficient slope, C A « -2.548 
Normal force coefficient slope, C n - -.7425 

n 

pitching force coefficient slope C m * -.31 

n 


Equations of Motion 

The EMS simulation equations were written only for the purpose of 
lunar return entry training for altitudes from 300,000 to 100,000 feet 
and therefore could be simplified to allow them to be programed on the 
three AD-80 analog computers which were available for the job. Although 
the two DDF24 digital computers were in use at the Flight Acceleration 
Facility, they were committed to control and monitoring functions and 
could not be used in the simulation of the entry vehicle. Even though 
some rather severe simplifications were made to the equations to fit them 
to the available computers, events proved that the equations were more 
than adequate for the objectives of the program. 

Coordinate Systems .- Two coordinate systems were used in the 
reentry equations. The (l) inertial coordinate system, and (2) CSM 
or body coordinate system are defined in the following two paragraphs 
ana are diagramed in figures 4 and 3 respectively. 

The inertial coordinate system is an orthogonal axes system whose 
origin is coincident with the surface of the earth. The inertial X T axis 
is parallel to the earth's surface and is positive pointing west. The 
Yj axis is perpendicular to the X^. axis, parallel to the earth's surface, 

and is positive pointing north. The Z^. axis is positive pointing toward 

the earth's center and completes the right-hsnded triad. 

The body cooYdinate system is an orthogonal system which is fixed 
to the vehicle and whose origin is on the command module longitudinal 
axis. The axis is coincident with the vehicle longitudinal sxis and 
is positive pointing through the apex of the command module. The Y^ axis 
is perpendicular to the X axis snd is positive pointing through the pilot's 
right arm. The axis is positive pointing through the vehicle center of 
gravity and completes the right-handed triad. 
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Mathematica l Model.- The math model (reference 2) used for this 
simulation was developed by Control Requirements Branch, Guidance end 
Control Division. The equations which are presented in block diagram 
form in figure 5 are written in six degrees of freedom and were simplified 
because of (l) the limitation of analog equipment, (2) limited time 
available to mechanize the equations, and (3) limited objectives (simula- 
tion was to be used only for astronaut crew training). However, the 
equations were sufficiently detailed to correctly calculate the vehicle 
dynamics and accelerations necessary to simulate lunar entries from .05 g 
(300,000 ft altitude) to 100,000 ft altitude. In addition to the normal 
spacecraft equations of motion, the EMS range-to-go equations (reference 3 ) 
were mechanized and solved on the computers. 

The principal simplifications and assumptions made in writing the 
equations were as follows: 

a. The earth model is a non-rotating cylinder with a constant 
gravitational force. 

b. Spacecraft position calculations are omitted except for altitude. 

c. Cross -coupling Inertia and products of inertia terms are omitted. 

d. RCS jet interaction terms have been omitted. 

e. Center of gravity offset in the Yg axis has been neglected. 

f. Aerodynamics are functions of total angle of attack only and 
are considered linear functions of the cosine of the angle of attack. 

g. Mass and inertia characteristics are constant. 

Special mention should be made of the first assumption. A centrifugal 
force term is subtracted from the gravitational force and this difference 
is multiplied by a gravity gain factor (K-, in the equation). This 
factor was adjusted experimentally so that the duration of the runs and 
the maximum acceleration were approximately the same as those obtained 
in more sophisticated simulations. 

Inputs to the equations of motion were the commanded rotational 
rates from the pilot's three axes RHC (rotational hand controller in 
the cockpit). 

The analog outputs (solutions of the equations) were recorded on 
two eight-channel strip chart recorders. In addition, three of these, 
the spacecraft attitude, EMS roll stability attitude, and the vehicle 
applied accelerations, were sent to the Honeywell digital computer. 
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Analog Mechanization 


The equations of motion mechanized on the analog computers are 
shown In a detailed block diagram in figure 5 . A majority of the 
mechanization was of standard configuration and no attempt will be made 
to explain the entire program. However, some explanation will be given 
on three of the circuits: 


a. Generation of sing^, cos 0 , V , and cos tj . - The technique 

used to compute the trigonometric functions of the aerodynamic angles 

and the total velocity was based on a paper by Robert M. Turner (reference k ) . 

It is an implicit or bootstrap m ethod whi ch solves an equation of the 

type R « XcosG + YsinG where r » + y£ end © « tan"* 5 - Y/X rather than 

the more straightforward computations R « ^x 2 + Y 2 , sinO « T sssspsasss as^ and 
X 


cobG 


m ttn— 

VX* 
‘cults an 


DM 

Y* 


The implicit technique eliminates the square root 


circuits and places the inverse trigonometric functions inside e closed 
loop. Experience has shown that this technique is more accurate and 
more stable than the more conventional technique when the total angle 
of attack approaches zero. When this happens, both X end Y approach 
zero and the conventional circuits run into trouble by attempting to 
divide one small quantity by another small quantity. Circuit diagrams 
are shown in figure 6. 


b. Generation of dynamic pressure, q . - The atmospheric density, p , 
varies by a factor of nearly 10^ for the altitude range encountered 
during an entry. If £ were generated directly then at an altitude above, 
say 200,000 feet, the signal-to-noise ratio would be so low that accurate 
generation of dynamic pressure would be impossible. This difficulty was * 
avoided by generating p u, squaring it, multiplying this variable by 

velocity, and finally squaring p i V to obtain dynamic pressure. To use 
the full range of the function generator, the function (2 pi - f i ) 

versus (h - 200,000) was programed. Details of mechanization are given 
in figure J, 


c. Mode control of the analog integrators . - For checkout and opera- 
tional purposes, the three analog computers were slaved together with 
one computer acting as the master. Under these conditions the operational 
modes, initial conditions (ic), hold (h) and operate (CP), of the slave 
computers are synchronized to those of the master computer. For ease end 
convenience of checkout and opera-tion, the mode control signals were 
generated at two different sources, the master analog computer and the 
DDP2^ digital control computer with the control signals being gated to the 
analog computers by a switch on the master analog console logic board. The 
operational philosophy employed was to use the master analog control signals 
during open loop testing (unit checkout) and to use the digital control 
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signals during the enalog-digltal-ccntrifuge closed loop runs. The logic 
circuitry required for the options} control of the analog integrators is 
shown in figure 8. 

Hun Termination .- The aforementioned mode control signals were also 
ured to terminate the data runs, that is, put the integrators into hold 
mode. There were two, normal and abnormal, types of terminate criteria 
to stop a close loop run. The normal run termination took place when 
the vehicle altitude reached 100,000 feet. Abnormal termination occurred 
when an analog amplifier overload lasted longer than one second, or a 
hardware failure other than the analog computers occurred, such as digital 
hardware failure, centrifuge overspeed, major power failure, or if the 
medical officer or pilot determined the run ought to be stopped. 

The altitude and analog amplifier overload terminate signals were 
generated on the analog and sent to the digital control computer where 
they were combined with the other terminate signals, and the digital 
control computer in turn sent a hold signal back to the analog integrators. 
At the same time, the digital control computer also sent control signals 
to the centrifuge to decrement the centrifuge arm speed to a level of 
1G at a rate of l/l6 G per second. The circuits used to generate these 
termination signals are shown in figure 8. 


RKSULTS AND DISCUSSION 

The EMS simulation was used to train the Apollo 8 prime crew 
(Borman, Lovell, Anders) and backup crew (Armstrong, Aldrin, Haise) 
for two days - 1 day in shirt sleeves and 1 day in vented suits, the 
Apollo 10 command pilot (John Young) for one day; and one additional day 
using two pilots (R. Lindsay and M. Lake) from Flight Crew Support Division 
was used to determine if entries could be made using a limited amount of 
visual cues (tha't is, no EMS instrumentation). It was demonstrated that 
safe lunar return entries could be made using only a G meter and a roll 
attitude indicator as visual cues. 

The initial conditions for the various entry profiles that were flown 
to take data are given in table 1 with all runs beginning at approximately 
300,000 feet with the lift vector up. During each run the command module 
stayed in free flight from 300,000 fteet to .05 g's (^290,000 ft) at which 
time the pilot took over the controls and flew to a precalculated EMS 
target at 100,000 feet altitude. Two techniques (l) EMS ranging and 
(2) constant g profiles were used to fly to the precalculated target. The 
success of the run was determined by how close to zero the RTG counter 
reed at the end of the run, l*e., when altitude 100,000 feet. An 
example of the analog recorded results from a constant 3 g profile nan is 
shown in figure 
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During the sctusl data runs, the simulation flew wi J, i no analog 
or digital computer failures. Troubles that did occur during the data 
runs occurred in the cockpit instrumentation, (l) EMS RTG counter, 

(2) the G-meter, and ( 3 ) roll stability indicator (RSl) meter. These 
Instrument troubles were minor and did not hamper the operation of the 
simulation in taking data. The troubles were intermittent and were of 
the type where the instrument needle or indicator did not fellow the 
ana log -digital signal tha'c was sent to them. 


Digital Verification of the EMS Math Model 

To validate the EMS equations of motion, a digital simulation using 
a fourth order Runge-Kutta integration routine with a .15 sec. step size 
was run and the results compared with a supercircular entry run from the 
CSM Guidance and Control Verification Simulation (CSM G&CVS). The EMS 
equations of motion were considered verified when the plots of altitude, 
total velocity, and dynamic pressure versus time followed the same general 
pattern of those from the CSM G&CVS. 

Following this, additional digital runs were made to obtain the proper 
gravity coefficient gain factor, Kp (required in the EMS equations to 
compensate for the non-rotating earth). The gain factor was determined 
by the agreement of the EMS plots of altitude, dynamic pressure, and 
velocity versus time with similar plots from the CSM G&CVS simulation. 

A Kq = 1.06 was required and the compared plots are shown in figure 10. 

A further check on the equations of motion was performed when the 
digital program was used to generate range -to-go (RTG) numbers for a 
constant G range-to-go chart. To obtain the RTG numbers, additional 
guidance equations (reference 3 ) were added to the digital program to 
make the vehicle automatically fly constant G entry profiles. These 
additional equations shorn in figure 11 were used to obtain range-to-go 
numbers for constant 1, 3 , 4, and 5 G f s for the four cases used to take 
close loop data. These RTG numbers were charted for the various constant G 
runs and then used to initialize the EMS RTG counter prior to each data 
run. The EMS range-to-go results compared with MPAD's constant G entry 
RTG numbers within 0-100 n.m. for 800-1600 n.m. runs. Comparison of the 
two sets of RTG numbers were tabulated and are shown in table 2. 


Dynamic Check of the Analog Mechanization 

Three types of runs using lunar return entry initial conditions 
were used to dynamically verify the analog mechanization of the EMS 
equations. The success criteria for the analog mechanization was to 
compare the analog solutions of the three runs with digital solutions 
of the identical runs. The three types of runs were: 

1. No aerodynamics (zero dynamic pressure). 
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2. Constant bank angle (£ * -90°) with aerodynamics. 

3. Continuous 20°/s rolling entry with aerodynamics. 

The comparison of the digital and analog dynamic solutions of these runs 
are shown in figures 12, 13, and 14. 

Initially during the checkout period a considerable amount of trouble 
was encountered in obtaining satisfactory analog solutions for these runs, 
that is. to match the digital solutions. The trouble was discovered to be 
an offset in 9, rate of eh n nge of the Euler pitch angle which resulted in 
a drift in the vertical velocity. The offset was constant for a particular 
day but would change from day to day. Initially (when it was first discovered), 
the 5 offset was 0 . 7 ? mil lira dian/see but this was due to a bad multiplier in 
the circuit forming 9. This was replaced and the offset dropped to 
0.10 - 0.30 milliradisns/sec depending on the day. To compensate for this 
offset, a bias voltage was '-dded to the Q circuit ry which temporarily 
eliminated the drift. However, long term changes in the multiplier offset, 
which were within normal drift tolerances, required that a careful and 
almost continued check be made on the sensitive circuit. This problem 
is discussed in more detail in reference 5 and in the appendix of this 
paper. 


CONCLUDING REMARKS 

A simplified set of equations was used in a closed loop entry 
simulation to train the Apollo 8 prime and backup crews. The equations 
of motion were solved on analog computers, the outputs of which were 
used as drive signals for the Flight Acceleration Facility centrifuge. 

The digital solutions of the simplified equations and the actual piloted 
data, runs verified -the equations to be adequate for the purpose of training 
pilots for lunar entries. The digital verification results of the EMS 
equations agreed with results obtained from the CSM G&CVS simulation and 
from an MPAD digital entry simulation. In addition to training the 
Apollo 8 crews, the simulation was used to demonstrate that’ safe entries 
could be made from lunar return velocities using only a meter for an 
indication of roll angle and an accelerometer as instruments. 

During the actual data runs the simulation operated with no failures 
or major problems with the analog or digital computers. However, during 
the checkout phase a ^ sensitivity problem was encountered in solving the 
equations of motion. This resulted in a drift in vertical velocity which 
was eliminated only by a careful and almost continued check on the sensitive 
circuit. The equations have been modified in an attempt to reduce or 
eliminate the sensitivity and it is recommended that the modified equations 
be programed to determine if the sensitivity has, indeed, been eliminated. 
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APPENDIX 


Modifications to Equations of Motion for Entry Simulation 


The equations of motion used in the entry simulations seemed to be 
unusually sensitive to small errors in the computation of 9, the derivative 
of the pitch Euler angle. The reason for this can be seen by the following 
analysis. 


Starting with the equations as given in figure 5 end assuming that the 
yaw angle, yaw rate, lateral velocity, rolling velocity, and roll angle 
and dynamic pressure are all zero, the .equations of motion can be written 
as: 


u = -g e sin9 


w = g cos9 
°e 

e = £ 


h = u sin 9 - w cos 9 

e e = g - v t 2 /r 

where u and w are the x~ and y-body components of the velocity, 9 is the 
pitch angle, h is the altitude, £ is an assumed error in the computation 
of 9, V,j) is the total velocity, g the gravitational acceleration and H 
the radius of the earth. 


Taking the derivative of fi with respect to time, 

m 

h = u sin 9 - y cos 9 


+ 6(u cos 9 + w sin 9) 


which simplifies to 

h = -g^ + 9(u cos 9 + w sin 9). 

But u cos 9 + w sin 9 is approximately equal to the total velocity, 
and thus 1 

fi = -g e + e v T . 

• 

In the entry simulation, the static offset in 9 was about 0.75 milli- 
radians/second and the velocity was about 36,000 feet/second. The 
resulting error in vertical acceleration was 27 ft/sec^, the same order 
of magnitude as the gravitational acceleration. 
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If the equations are written so that the translational acceleration-to- 
velocity integrations are done in a non-rotating coordinate system, this 
problem would be eliminated . A suitable set of equations follows without 
derivation. The assumptions inherent in these modified equations are 
basically the same as for the original equations. In addition, simplified 
equations were used to compute the angles of attack and sideslip and it was 
assumped that the lateral acceleration was negligible. This was done to 
reduce the number of coordinate transformations, No changes were made in 
the equations for the Entry Monitor System or the control system. 


i 


I 


I. Aerodynamic -force equetiona 
A = -C A q S/m 


» 0 


\ “ - C N* S /“ 


II. Angular accelerations 

* = c ?££ ?§a + h. 

P Y D I + I 


C + c i ass + 

C m D C N I + I 

y y 


f - C^j qsr/l z + Kj/I 


III. Euler rates - pitch, yaw, roll . sequence 

V * q sin 0 + r cos 0 
© ** (q cos 0 - r sin 0)/cos^ 

0 * p - © sin 




•9fam 
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IV. Aerodynamic accelerations, inertial axes (assume A ^ * 0) 

A * A cos ty cos © + A (cos 0 sin © + sin 0 sin^/ cos ©) 
x i *b h 


A * A sinl/^ - A sin 0 cos ^ 
y i *b Z b 

A E -A sin 9 cos I/' + A (cos 0 cos © - sin 0 sin^ sin 0) 


V, Total inertial accelerations 


V ** A 
x x, 


V * A 

y y t 


•\ b + (¥ - *) K i 

^ + (t- • j h 


VI. Miscellaneous velocities and angles 


2 + v M 


V H- (V * + V y> 


2 + V 2 }^ 


V T = (V H + V> ; 


= flight path angle 


gnt 


u 


VII. 


« heading angle 
- tan' 1 (-V /-V x ) 

* eln" 1 ( -Vy/V H ) 

* angle of attack 

* tan - ^ (w/|u| ) 

2 - |(© ♦ y)coB0 ♦ {v - t/ h ) R in 0 

^ «= angle of sideslip 

“ sin' 1 (v/V T ) 

Z -Re *'f)oin 0 - (V- )cos0| 


u 


h 

X 2 

X 3 


V 

X 

V 

s 

m i • 

m 2 

"3 

0 

V 

y 

V 


n i 

C\J 

a 

n 3 


v l 

z 


Aerodynamic coefficients snd dynamic pressure 


C, c axial force coefficient 

A 


C A + C 1 * 

M o 1 


* normal force coefficient 


C »0 * V* 


C Y = lateral force coefficient 


c„- p 




C ■ pitching moment coefficient 
m 


C ♦ C 
nu 








IX. Constants 


l^,» lyl ij, t ^ ^ 

C . t C. j C M , C-j , Cy t c , c 
A 0 1 N o N « V m O m <X 

C «y V M J’ V 
Function of altitude 

* 

p ■ mass density of air 


The modified equations have been programed and checked out for the 
three ADSO analog computers used in the original simulations and a 
digital check simulation has been obtained. The results of these 
simulations have been compared with results obtained from more elaborate 
simulations such as the Command Module G&C Verification Simulation and 
generally good agreement has been obtained. Figure 15 is a comparison 
of results for a rolling entry obtained from the simplified analog 
simulation and results obtained from the Command Module G&C Verification 
Simulation study. As an inspection of the figure shows f there is good 
agreement between the two simulations. The differences in the dynamic 
pressure probably arise from differences in the roll rate or other minor 
differences in the control system. Although not shown on the figure, the 
short period dynamic characteristics are also in good agreement. The 
same control system was used in the modified equations as in the original 
equation and hence the control characteristics which were reported to 
be good in the original simulation should be equally good with the 
revised equations. 

The modification to the equations have removed the undue sensitivity 
encountered in the original equations and no new undesirable characteristics 
have been experienced. These facts coupled with the good agreement with 
the hybrid entry simulation make it advantageous to use the modified equa- 
tions in future entry simulations where there are severe restrictions on 
the amount of computing equipment. 


C n - yawing moment coefficient 



q * dynamic pressure 

-J/OVt 2 

VIII. Direction cosines 

1^ * cos 9 cos p 

1 2 » sin \y 

1^ * -sinOcos 

m 1 * sin 9 sin 0 - cos 9 sin V cos 0 

Bg « eosi^ cos 0 

*■ cos 9 sin 0 + sin 9 sin*/' cos 0 

*= sin 9 'bos 0 + cos 9 sin *y sin 0 

n^ * -cos {/ sin 0 

* cos 9 cos 0 - sin 9 sin sy sin 0 

NOTE: Only six of the nine direction cosines are needed. The 

direction cosines m^, nig, and sre not required. 




TABLE 1 INITIAL CONDITIONS TOR CLOSE LOOP DATA RUNS 
















































Flight path a . .lei 

aT 400000 FEET 
(OEG) 


CONSTANT g 
PROFILES 
(G'S) 


RANGE TO GO FROM .05G'S. 
^ 290000 FEET ALTITUDE,, 
TO 100000 FEET ALTITUDE 

(NAUTICAL MILES) 



TABLE 2. EMS DIGITAL RTG VALUES COMPARED WITH THOSE 
FROM AN MPAD (MISSION PLANNING AND ANALYSIS DIVISION) 
DIGITAL ENTRY SIMULATION FOR ENTRIES FROM .OSG'S TO 
100000 FEET ALTITUDE. 


































IGURE 1 BLOCK DIAGRAM OF TOTAL SYSTEM 









FIGURE 2 BLOCK DIAGRAM OF RCS RATE COWAND CONTROL SYSTEM. ROLL 
CHANNEL SHOWN. PITCH AND YAW CHANNELS ARE SIMILAR. 
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FIGURE 4 DIAGRAM SHOWING ORIENTATION OF INERTIAL COORDINATE STSTEM 






FIGURE! 5 BLOCK DIAGRAM Of If L> MAM MODE. 

















FIGURE 6 ANALOG CIRCUIT USED TO GENERATE SlNE(C)p ), COS I"* 1 " V../ ,^uo I NE(T|) ,V T 

























9 - continue! 
















1. CONSTANT -G guidance equations and logic 


2. ROLL RATE COWAND LOGIC 

3. RANGE TO GO EQUATIONS 


1 


VT m 

rAklCT ImT a 

Ax f A ’ ^ 

LUN j 1 AN 1 -U 

GUIDANCE equations 

1 AT 1 r 


- PD - 

A» 4U LvAj 1 V- 



WHERE : 


jiL. 


iLUiL 


RANGE TO GO] 
EQUATIONS 



V T -TOTAL VEHICLE ACCELERATION 
Ax,Az -APPLIED ACCELERATIONS 
fl -ALTITUDE RATE 
DD -DESIRED DRAG 

Qvc -LIFT VECTOR COWAND ATTITUDE 
Pc -ROLL RATE COWAND 
Vi (»c)*Vxb AT .C6G'S 



FIGURE 11 ADDITIONAL EQUATIONS USED TO GENERATE EKS RANGE T0 GO 
VALUES FOR CONSTANT-G ENTRY PROFILES. 
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figure lU - Dynaal: eJwck ca*e 3, continuous 20 °/s rolling entry vlth •erodyna»lc« 
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